Critical phenomena in atmospheric precipitation 
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Critical phenomena near continuous phase transitions are typically observed on the scale of wave- 
lengths of visible light Here we report similar phenomena for atmospheric precipitation on scales 
of tens of kilometers. Our observations have important implications not only for meteorology but 
also for the interpretation of self-organized criticality (SOC) in terms of absorbing-state phase transi- 
tions, where feedback mechanisms between order- and tuning-parameter lead to criticality. 0] While 
numerically the corresponding phase transitions have been studied, 0,13 we characterise for the first 
time a physical system believed to display SOCQ in terms of its underlying phase transition. In 
meteorology the term quasi-equilibrium (QE)Q refers to a state towards which the atmosphere is 
driven by slow large-scale processes and rapid convective buoyancy release. We present evidence 
here that QE, postulated two decades earlier than SOC^, is associated with the critical point of a 
continuous phase transition and is thus an instance of SOC. 

Journal reference: Nature Physics 2, 393 - 396 (2006). doi:10.1038/nphys314 
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Self-organized criticality has been proposed as an ex- 
planation for scale-free behaviour in many different phys- 
ical systems 0. In most of these, however, it is im- 
possible to measure standard observables for critical 
phenomena, such as order parameters, tuning parame- 
ters, or susceptibilities. Consequently, despite theoretical 
advances^, IE SOC has only loosely been connected 
to the broader field of critical phenomena. The present 
study helps position SOC as a sub-branch of critical phe- 
nomena by examining a system where the identification 
and measurement of standard observables is feasible. 

At short time scales the majority of tropical rainfall oc- 
curs in intense rain events that exceed the climatological 
mean rate by an order of magnitude or more. Precipita- 
tion has been found to be sensitive to variations in water 
vapour along the vertical on large space and time scales 
both in observations 9, 10] and in models . [ill ll^ fl^ This 
is due to the effect of water vapour on the buoyancy of 
cloud plumes as they entrain surrounding air by turbu- 
lent mixing. We conjecture that the transition to intense 
convection, accompanying the onset of intense precipita- 
tion, shows signs of a continuous phase transition. The 
water vapour, w, plays the role of a tuning parameter 
and the precipitation rate, P(w), is the order parameter 
Note that such a large-scale continuous phase transition 
involving the flow regime of the convecting fluid is en- 
tirely different from the well-known discontinuous phase 
transition of condensation at the droplet scale. We ana- 
lyzed satellite microwave retrievals of rainfall, P, water 
vapour, w, cloud liquid water and sea surface tempera- 
ture (SST) from the Tropical Rainfall Measuring Mission 



from 2000 to 2005. Observations from the western Pacific 
provided initial support for our conjecture: a power-law 
pick up of the order parameter above a critical value of 
the tuning parameter, Wc, was observed. We preceded to 
test whether other observables also behaved as predicted 
by the theory of phase transitions. 

As motivation for our conjecture consider a generic 
lattice-based model which exhibits a continuous phase 
transition. Particle-conserving rules defining the model 
ascribe a number of particles to every lattice site, and 
demand hopping of particles to nearest-neighbour sites 
when a local density threshold is exceeded. The global 
effect of these rules is a phase transition at a critical 
value of the global particle density between a quiescent 
phase (where the system eventually settles into a stable 
configuration) and an active phase (where stable config- 
urations are inaccessible). The timing parameter is the 
particle density and the order parameter is identified as 
the density of active sites [14j. 

SOC can be described in terms of such absorbing-state 
phase transitions. 0, Q Here a coupling between order 
parameter and tuning parameter is introduced by open- 
ing the boundaries and adding a slow drive: whenever 
activity ceases, a new particle is added to the system, 
i.e., an increase in the tuning parameter. Large activ- 
ity on the other hand leads to dissipation (particle loss) 
at the boundaries, i.e., a reduction of the tuning param- 
eter. Such open, slowly driven systems organise them- 
selves to the critical point of the corresponding (closed 
boundaries, no drive) absorbing state phase transition. 
The critical behaviour as derived from finite-size scaling 
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analyses is the same in both cases |3, 14 , although the 
reason for this universality is not fully understood. [isj 
The scale-free avalanche size distributions in SOC mod- 
els result from the proximity of the system to a critical 
point. 

From the meteorological perspective, a related moti- 
vation for our conjecture arises. Atmospheric convec- 
tion has long been viewed similarly in terms of a slow 
drive (surface heating and evaporation) and fast dissipa- 
tion (of buoyancy and rainwater) in precipitating con- 
vection. Surface heating and evaporation drive turbu- 
lent mixing that maintains a moist atmospheric bound- 
ary layer. Combined with radiative cooling, conditional 
instability is created — while sub-saturated air remains 
stable, saturated condensing plumes can rise through the 
full depth of the tropical troposphere. The fast dissipa- 
tion by moist convection prevents the troposphere from 
deviating strongly from marginal stability, [ig Although 
observational tests of this approximate QE state of the 
tropical troposphere have limited precision, it forms the 
basis of most convective parameterizations in lar ge s cale 
models ^3 ^^'^ much tropical dynamical theorv.|l^ Hgf 
Taking large-scale flows into account modifies the pro- 
cess in space and time but does not change it fundamen- 
tally. This perspective suggests that a critical point in 
the water vapour would act as an attractor. Indeed this 
is basically the convective QE postulate.]^ 

The critical value Wc depends, e.g., on atmospheric 
temperature, but for present purposes this translates well 
enough into a critical amount of water vapour for a given 
climatic region. Regions here are defined by longitude 
ranges given in the caption of Fig. 1 corresponding to 
major ocean basins, for oceanic grid-points within 20S- 
20N. Data are collected at 0.25 degree latitude-longitude 
resolution. The observable w captures vertically inte- 
grated, or column, water vapour. It is given as a volume 
per area in units of mm. 

In Fig. 1 we show as a function of the tuning parameter 
w the average value of the order parameter (P) (w) and 
the susceptibility of the system, represented by the order 
parameter variance, ap{w), discussed following Eq. |(2Jl. 
The ensemble size for the average ranges from a few thou- 
sand at extremes to 10^ at typical ly-values. Above Wc, 
the order parameter is well approximated by the stan- 
dard form I 



(1) 



where a is a system-dependent constant and /? is a univer- 
sal exponent. The deviations from power-law behaviour 
below Wc in the main graph of Fig. 1 are typical of critical 
systems of finite size.|2l| 

The critical value Wc is non-universal and changes with 
regional climatic conditions, as does the amplitude a. To 
test the degree to which curves from different regions 
i collapse, we re-scaled the w-values in Fig. 1 by fac- 
tors reflecting the non- universality of Wc and (P) (w) 
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FIG. 1: Order parameter and susceptibility. The 

main figure shows the collapsed (see text) precipitation rates 
(P) (it)) and their variances ap{w) for the tropical Eastern 
(red) and Western (green) Pacific as well as a power-law fit 
above the critical point (solid line). The inset displays on 
double-logarithmic scales the precipitation rate as a function 
of reduced water vapour (see text) for Western Pacific (green, 
120E to 170W), Eastern Pacific (red, 170W to 70W), Atlantic 
(blue, 70W to 20E), and Indian Ocean (pink, 30E to 120E). 
Data are shifted by a small arbitrary factor for visual ease. 
The straight lines are to guide the eye. They all have slope 
0.215, fitting the data from all regions well. 



and (Jp{w) by fp and /^2, respectively (setting Western 
Pacific factors to one). For visual clarity, the data col- 
lapse in Fig. 1 is shown only for the Eastern and West- 
ern Pacific — climatically very different regions. Similar 
agreement occurs for other regions (steps in the rescaling 
and figures for all regions are provided in the Supplemen- 
tary Information) . The exponent (3 seems to be universal 
and independent of the climatic region. In the inset to 
Fig. 1 we show the average precipitation as a function 
of the reduced water vapour Aw = (w — Wc)/wc in a 
double- logarithmic plot. Importantly, power laws fitted 
to these distributions all have the same exponent (slope) 
to within ±0.02. The data points in Fig. 1 represent the 
entire observational period, including all observed SSTs. 
Conditioning averages by SST ranges yields similar re- 
sults (see Fig. 3 and Supplementary Information), reduc- 
ing the subcritical part of the curves slightly. 

We define the susceptibility x(w^; L) via the variance of 
the order parameter Cp: 



x{w;L)=L''aUw;L), 



(2) 



where d denotes the dimensionality of the system and L 
the spatial resolution. Fig. 1 shows a suggestive increase 
in (Tp near Wc, and indicates that standard methods for 
critical phenomena can sensibly be applied. 

Next we test for finite-size scaling. Because our system 
size cannot be changed, we identify the spatial data res- 
olution L as the relevant length scale. Changing L has 
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FIG. 2: Finite-size scaling. The variance of the order pa- 
rameter a%{w) as a function of w, rescaled with L° *^ for 
system sizes 0.25°, 0.5°, 1°, and 2° in the Western Pacific. 
From u) « 57 mm, this produces a good collapse. The inset 
shows that away from the critical point, up to w ~ 40 mm a 
trivial rescaling with L'^^^ works adequately. This suggests 
that the non-trivial collapse is indeed a result of criticality. 
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FIG. 3: Residence times. The number of times N{'w) 
an atmospheric pixel of 0.25° x 0.25° was observed at water 
vapour w in the western Pacific, given a sea surface temper- 
ature within a 1°C bin at 30° C. The green and blue lines 
show residence time for all points and precipitating points, 
respectively. The red line shows the order-parameter pick-up 
(P) (w) for orientation (precipitation scale on the right). 



the effect of taking averages over different numbers of de- 
grees of freedom and allows one to investigate the degree 
of spatial correlation. The finite size scaling ansatz for 
the susceptibility is 

x{w;L)^L'y/''x{AwL'^l, (3) 

defining 7 and v as the standard critical exponents and 
the usual finite-size scaling function xi^)j constant for 
small arguments \x\ «C 1 and decaying as |a;|~^ for large 
arguments \x\ ^ ^ UM The variance ap{'w;L) is affected 
by uncertainties in w and Wc, making precise quantifica- 
tion of x(u'; L) difficult. We therefore do not estimate 7 
from the w-dependence of i), corresponding to large 
arguments \x\ in Eq. ^ but effectively fix Aw = — x 
and obtain 7/1/ from the L-dependence of the maximum 
susceptibility x'^'^^iL). 

The variance of the average (P) {w; L) over L'^ inde- 
pendent degrees of freedom decreases as crp(w; L) oc L~'^. 
In a critical system, however, the diverging bulk correla- 
tion length ^ oc {Aw)~^ ^ L (small argument in Eq. 
prohibits the assumption of independence. In this case 
Eq. lO with Eq. |2Jl yields 

cr|,'"^'^(i) OC L^/'^^'^. (4) 

Coarsening the spatial resolution of the data, we find in 
Fig. 2 that ap'°''''(L) scales roughly as L^^, with A = 
0.46(4). This suggests the exponent ratio ^/v ~ 1.54(4). 
At criticality, the spatial decay of correlations between 
order parameter fiuctuations becomes scale- free, [ioj This 
is equivalent to a non-trivial power-law dependence of 
the order-parameter variance on L (see Supplementary 
Information for details and conditions). Hence, Fig. 2 



indicates a scale-free correlation function of fluctuations 
in the rain rate in the range of 25 km to 200 km. This sug- 
gests that the meteorological features known as mesoscale 
convective systems ji^l are long-range correlation struc- 
tures akin to critical clusters. '24'| Synoptic inspection in- 
dicates that the high rain rate phase and critical region 
of Fig. 1 come substantially from points within such com- 
plexes (examples are provided in the Supplementary In- 
formation) . The question of self-organisation towards the 
critical point of the transition is addressed by display- 
ing the residence times of the system in Fig. 3. This is 
the number of observations in the 5-year period where 
the system was found at a given level of water vapour. 
A slowly driven system would be expected to spend a 
significant amount of time in the low-w phase because 
when it fluctuates into this phase e.g. due to some large- 
scale event, it takes a long time to recover. Therefore 
the distribution decreases slowly towards low values of 
w. The fast dissipation mechanism, on the other hand, 
ensures that the system leaves the high-w regime rela- 
tively quickly when it fluctuates into it. Consequently 
the distribution decreases rapidly towards large values of 
w. For the properties of rainfall, the part of the distribu- 
tion in Fig. 3 comprised only of observations with rainfall 
is of interest, seen as the blue line in Fig. 3. We note that 
the system is most likely to be found near the beginning 
of the intense precipitation regime. Almost the entire 
weight of the distribution of rainy times is concentrated 
here. 

Meteorologically, these results suggest a means to re- 
define and extend convective QE, both empirically and 
theoretically. In its simplest application QE assumes that 
the relationship among atmospheric column thermody- 



4 



namic variables is pinned close to the point where deep 
convection and precipitation set in. Fig. 3 shows this 
to be a reasonable first approximation , but it also im- 
plies associated critical phenomena. A loss term (P) {w) 
of the form of Eq. implies the absence of any well- 
defined convective time scale. Scale-free distributions of 
event sizes Q and the spatial correlation behavior seen 
in Fig. 2 may result from this proximity to an apparent 
continuous phase transition. 

These findings beg for a simple SOC-type model 
of the atmospheric dynamics responsible for the crit- 
ical behaviour. While the physics must conform 
with recent cloud-resolving model analysis of mesoscale 
aggregation [ill E^l]. our results point to the key role of 
excitatory short-range interactions, essential for critical 
phenomena of the type seen here. This study advances 
our understanding of SOC as a critical phenomenon, 
identifying the underlying phase transition and associ- 
ated critical phenomena. Beyond scale-free event size dis- 
tributions it furnishes direct evidence, for the first time, 
for an underlying phase transition in a physical system. 

Methods — Data are from the TMI (TRMM mi- 
crowave imager), processed by Remote Sensing Systems 
(RSS). The spatial resolution reflects the footprint of the 
instrument. As with any satellite retrieval product, it is 
necessary to consider whether the algorithm assumptions 
could impact the results. The microwave retrieval algo- 
rithm is that used on Special Sensor Microwave Imager 
(SSM/I) data. The combination of four microwave 
channels permits independent retrieval of water vapour 
and condensed phase water (with SST and surface wind 
speed), while an empirical relation is used to partition 
cloud water and rain. Column water vapour validates 
well against in situ sounding data, which also show that 
daily variations are largely associated with the lower tro- 
posphere above the atmospheric boundary layer. Val- 
idation of TMI rain rate against space-borne precipita- 
tion radar (PR) at sub-daily time scales in the tropical 
Western Pacific [l^ show TMI overestimating rain rate 
but with an approximately linear relationship to PR. We 
have performed a number of checks to verify that results 
are not substantially impacted by a high rain rate cutoff 
in the algorithm (25 mm/h), including comparison to re- 
gions where cutoff occurences are very low, such as the 
eastern Pacific (Fig. 1). The clearest check is that the 
essential features are identical for the cloud liquid water, 
whose measurement cutoff of 2.5 mm is never reached (see 
Supplementary Information). SST data here are averages 
over non-flagged neighbors in space and time, since SST 
is not retrieved at high rain rates. 

The critical value Wc is determined by an iterative pro- 
cedure with an initial guess, followed by a fit to Eq. 
above Wc- Error bars in Fig. 3 are standard errors of 
the variance crp(w; L), determined via the zeroth, second 
and fourth moments of P{w). Individual measurements 
of P(w) are considered independent, which holds well be- 



tween satellite overpasses, though not within individual 
tracks. 
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